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ABSTRACT 

(-H \ Wc study how to recover the fuU 3D clustering information of P(fc, z), including 

^ J redshift space distortions (RSD), from 2D tomography using the angular auto and 

cross spectra of different redshift bins C(,{z, z'). We focus on quasilinear scales where 
the minimum scale Xmin or corresponding maximum wavenumber kmax = '^'^ I ^min is 
targeted to be between kmax = {0.05 — 0.2}/iMpc~^. For spectroscopic surveys, we 
find that we can recover the full 3D clustering information when the redshift bin width 
IS.Z used in the 2D tomography is similar to the targeted minimum scale, i.e. Az ~ 
{0.6 — 0.8} XminH (z) / c which corresponds to Az ~ 0.01 — 0.05 for z < 1. This value 
of Az is optimal in the sense that larger values of Az lose information, while smaller 
values violate our minimum scale requirement. For a narrow-band photometric survey, 
with photo-z error ct^ = 0.004, we find almost identical results to the spectroscopic 
survey because the photo-z error is smaller than the optimal bin width az < Az. For a 
typical broad-band photometric survey with = 0.1, we have that az > Az and most 
radial information is intrinsically lost. The remaining information can be recovered 
from the 2D tomography if we use Az ~ 2az- While 3D and 2D analysis are shown 
here to be equivalent, the advantage of using angular positions and redshifts is that we 
do not need a fiducial cosmology to convert to 3D coordinates. This avoids assumptions 
and marginalization over the fiducial model. In addition, it becomes straight forward 
to combine RSD, clustering and weak lensing in 2D space. 
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1 INTRODUCTION 

In recent years, galaxy redshift surveys have provided new 
information about the cosmological model of our Universe, 
in pace with precision cosmology from other probes like 
CMB and type la Supernovae. We are now entering excit- 
ing times for cosmology, when surveys will go deeper and 
wider with increasing number of galaxy positions in each 
catalogue. With deep surveys we can use weak lensing infor- 
mation to improve constraints on cosmological parameters 
and also trace directly the dark matter distribution at large 
scales. Theoretical analysis of weak lensing (WL) is usu- 
ally made through a 2D (angular) analysis of the measured 
galaxy shear maps. Future surveys will have less shot noise, 
allowing for more freedom in how we break the sample into 
multiple redshift shells, so that galaxy correlations can also 



be measured in and between shells. In doing angular corre- 
lations, we are projecting all the radial information within 
each redshift bin. But if we are able to use very thin radial 
shells, we can maybe recover the radial information using 
the angular cross-correlations between all the redshift bins 
(see Montanari & Durrer 2012 for a related idea). This is 
what we want to investigate in this paper. 

This goal is also connected to recent studies of galaxy 
surveys using a combination of redshift space distortions 
(RSD) and WL galaxy-shear and shear-shear correlations. 
These allow measurements of galaxy bias and the breaking 
of degeneracies between growth history and cosmic history, 
as has been recently proposed (Gaztaiiaga et al. 2012; Cai 
& Bernstein 2012). RSD are usually studied in 3D, which 
complicates a joint analysis with WL which is usually 2D 
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(see Kitching et al. 2011 for a comparative analysis with 
3D cosmic shear). If we could study RSD in 2D without loss 
of information, then it would be possible to do a joint anal- 
ysis of both probes using only angular correlations with the 
corresponding simplification in the covariance analysis. 

Observations directly probe redshifts and angular co- 
ordinates on the sky. Doing an angular analysis therefore 
does not require any prior knowledge of the cosmological 
model, while for doing 3D analysis we have to assume a 
fiducial cosmology to convert to comoving spatial coordi- 
nates. This then requires modelling the Alcock-Paczynski 
effect when fitting different models to our observables (Al- 
cock C, Paczynski B. 1979). As the transformation is red- 
shift dependent one has to make sure that this procedure 
is not biasing the parameter constraints. If the theoretical 
prediction for the correlations in angle and redshift can be 
calculated for each model, an angular analysis relating di- 
rectly to the observables is much more direct. 

The final goal of this paper is to analyze the bin op- 
timization that allows us to recover the 3D constraints on 
clustering using a 2D tomographic approach. We have stud- 
ied this in the framework of several idealized surveys: a spec- 
troscopic survey in a redshift range similar to SDSS redshift 
range; a survey with photometric redshifts from a camera 
with narrow-band filters like the camera that Physics of the 
accelerating Universe Survey (PAUj^i a^nd finally, a survey 
with redshifts obtained from photometry with broadband 
filters, in a redshift range similar to Dark Energy Survey 
(DESfl For the three surveys we have analyzed a bias fixed 
model, constraining r2,„. In addition, in the spectroscopic 
survey we have also studied the standard RSD constraints 
on the bias h and growth index 7. 

In section[2]we describe galaxy surveys, parameters con- 
sidered in the analysis and a description of the observables. 
In section [3] we show the constraints obtained in our anal- 
ysis for the different surveys described above. Finally, we 
summarize all the results in section [4] with the conclusions. 



Case n(r-) (h^Mpc^^) 



Low Shot-Noise 3.14 ■ 10"^ 

High Shot-Noise 6.89 ■ IQ-* 



Table 1. Comoving galaxy number densities sA, z = 0.55 assumed 
in this paper for the spectroscopic and narrow-band photometric 
surveys. Case 1 corresponds to a low shot noise level {nPgai ~ 2%, 
where Pgai is the monopole of the galaxy spectrum at z = 0.55 
and k = 0.1/iMpc~^) while Case 2 corresponds to a high shot 
noise level {nPgai = 10%) . 

Throughout the paper we used CAMB_source ^ (Challi- 
nor & Lewis 2011) to compute the matter 3D power spectra 
as well as the angular power spectra, including cross corre- 
lations between radial bins. 

2.1 Fiducial surveys and galaxy samples 

In this section we describe our fiducial surveys and galaxy 
samples. We characterize them by a redshift range, a given 
accuracy of redshift measurements, a galaxy redshift distri- 
bution and bias. 

In all cases we assume a full sky coverage. In ideal con- 
ditions this implies that the covariance matrix of observables 
such as Ci is diagonal in I (but notice that this assumption 
is not expected to change the conclusions of this paper). In 
all three surveys the overall redshift distribution of galaxies 
per deg^ is taken as, 

which is typical of a flux-limited sample with a magnitude 
cut at iAB < 24. In Eq. ([1} Ngai is a normalization related 
to the total number of galaxies per square degree under con- 
sideration. 



2 METHODOLOGY 

The goal of this paper is to show under which conditions, 
if any, one can recover the full 3D clustering information 
from a tomography study. By this we mean a combination 
of all the auto and cross angular spectra after the survey vol- 
ume has been divided in a set of consecutive redshift bins. 
The angular spectra within each bin will include information 
mainly from transverse modes, while cross correlations be- 
tween different shells accounts for radial modes with scales 
comparable to the bin separation. 

We investigate this idea in the context of a spectroscopic 
survey as well as two photometric surveys with different ac- 
curacies in the redshift determination. In what follows we 
describe these "typical" surveys, the assumed galaxy sam- 
ples, the observables considered and the figures of merit used 
to compare 3D and 2D tomography results. 

1 www.pausurvey 

2 www.darkenergysurvey.org 



2.1.1 Spectroscopic survey 

Our benchmark spectroscopic survey has radial positions 
given by true redshifts (i.e. CTz = in the formulation be- 
low) and a redshift range 0.45 < z < 0.65. Hence for the 2D 
tomography of this survey we use top hat binfl to compute 
angular power spectra. In Table [2] we show the different bin 
configurations considered, characterized by the number of 
bins in which we divide the survey volume and their width. 
Provided with the narrow redshift range we can assume that 
the bias does not evolve, hence we take b — 2 throughout. 

Lastly we discuss two cases for this survey, one where 
shot-noise is non-negligible and another where it is a sub- 
dominant source of error. These cases are detailed in Tabl^ 
and for the redshift range under consideration imply 9M and 
40M galaxies respectively (assuming full sky surveys). 

^ camb. info/sources 

^ To satisfy differentiability requirements at the edges we use in 
practice </i(^) oc exp [—((2 — z)/{Az/2))'^'^'\ , where z is the mean 
redshift of the bin and Az the full width. 
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Figure 1. Top panel shows the rcdshift distribution in the spec- 
troscopic and narrow band photo-z survey (violet). For the nar- 
row band case wo show how the true redshift distributions given 
by Eq. JTJ look like if we divide the volume in eight consecutive 
redshift bins. Bottom panel shows the same but for a broadband 
photometric survey divided in five bins. 



2.1.2 Narrow Band Photometric survey 

This case intends to be representative of a configuration such 
as the one proposed for the PAU survey where a set of nar- 
row band filters is expected to deliver "low-resolution" spec- 
tra in a redshift range actually broader than the one con- 
sidered here (Bemtez et al. 2009; Gaztanaga et al. 2012). 
Hence our narrow-band photo-z survey has accurate photo- 
metric redshifts of = 0.004, in the same redshift range of 
the spectroscopic case (0.45 < z < 0.65). The bias (fe = 2) 
and the shot-noise cases considered match those of Sec. 12.1.11 
(and are given in Table [T]). 

In turn the bin configurations assumed for the 2D to- 
mography are also the same as for the spectroscopic sur- 
vey given in Table [21 but with bin limits that now refer to 



Number of bins 


A2 


Ar {h 


-1 Mpc) 


1 


0.20 


4 


68 


4 


0.05 


113 


- 122 


8 


0.025 


56 


- 61 


16 


0.0125 


28 


- 31 


20 


0.010 


22 


- 25 



Table 2. Bin configurations used for the 2D tomography in the 
case of the spectroscopic and the narrow band photometric survey 
in a redshift range of 0.45 < z < 0.65. We show the number of 
radial bins and their range of widths in redshift and comoving 
distance. 

photometric redshifts. Thus the true redshift distribution of 
galaxies in each bin is no longer a top hat, but rather has 
a small overlap with the nearest neighbouring bins due to 
the photo-z error, as described in Eq. (|18|l below. In the top 
panel of Fig. [1] we show this effect for the particular case of 
8 bins. 



2.1.3 Broad Band Photometric survey 

On the other hand we consider a photometric survey that 
uses broad-band filters such as DES 0, Pan-Starrs [f| or the 
future imaging component of Euclid Q These surveys are 
expected to achieve photometric redshift estimates with ac- 
curacies ~ 5%/(l -I- z) (Banerji et al. 2008; Ross et al. 
2011). In what follows we do not consider a possible red- 
shift evolution of the photometric error but instead assume 
a conservative value of CTz =0.1. 

Typically optical photo-z surveys are fainter and sample 
a much larger number of galaxies than spectroscopic ones, 
hence we assume a broader redshift range, 0.4 < z < 1.4, 
and only a low shot-noise case as given in Table [T] For the 
redshift range assumed this implies ~ 150 x lO'' galaxies. 
Table [3] show the bin configurations we have considered for 
this case. While in the previous cases we have assumed the 
bias is constant with redshift (because of the narrow redshift 
range), for the broadband photometric survey we introduce 
an evolution following (Fry 1996), 

= l + (2) 

where 6* = 2 is the bias at = 1. In turn for the evolution 
of bias we have always assume the fiducial cosmology. 

2.2 Spatial (3D) power spectrum 

Since we are only interested in quasi-linear scales we as- 
sume the following simple model for the 3D galaxy power 
spectrum in redshift space, 

PAk,ix,z) = {h + fi,y D\z)Po{k)e-''"^"^''>^\ (3) 

5 www.darkenergysurvcy.org 
" pan-starrs.ifa.hawaii.edu 
^ www.euclid-imaging.net 
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Number of bins 


Az 


Ar {h 


-1 Mpc) 


4 


0.25 


398 


- 592 


5 


0.20 


315 


- 480 


6 


0.167 


260 


- 404 


7 


0.143 


221 


- 348 


8 


0.125 


193 


- 306 


9 


0.111 


171 


- 273 


10 


0.10 


153 


- 246 



Table 3. Bin configurations considered for a broadband photo- 
metric survey witliin a redshift range 0.4 < z < 1.4. We show the 
number of radial bins and their range of widths in redshift and 

comoving distance. 



where Po is the hnear spectrum a.t z = Q (properly normal- 
ized), D(z) is the linear growth factor and the remaining 
amplitude depends on the bias h{z) and the linear growth 
rate f{z) = d\n D / dlna. The Gaussian cut-ofT accounts for 
the fact that the radial information might be diluted due to 
photometric redshift errors 0. In Eq. ((3]) this redshift er- 
ror propagates to scales through ot{z) = ca^/ H (z). Notice 
that at depends also on the cosmic history. This should be 
taken into account when constraining relevant cosmological 
parameters (e.g. Sim). 

For a spatial analysis the measured 3D power spectrum 
depends on the cosmological model assumed to convert red- 
shift and angles to distances. Hence for every model to be 
tested against the data one must perform a new measure- 
ment. This process is very costly. Instead one can choose 
a reference cosmological model where the measurement is 
done once, and then transform the model prediction to this 
reference frame (Alcock C, Paczynski B. 1979). 

Let us call P°'"'{k,fi) the power spectrum measured 
in the reference cosmology and P'"°'^{k, fi) the model pre- 
diction at the point in cosmological parameter space being 
tested. The transformation of distances and angles from the 
cosmological model being tested {k, fi) to those in the refer- 
ence model (k, fi) is done through the scaling factors 



H{z) 



\Z) 



Cl_ = 



dmod 
A 



(^) 



dA^z) 



(4) 



as fcy = fc||/c|| and k± = k±/c±, where || indicates modes 
parallel to line of sight and _L perpendicular. The Hubble 
parameter and the angular diameter distances are given by 



H{z) = 100fe^n„(l-H2)3 + f7l3E(H-z)-3(l + "') (5) 

c dz' 
Jo -tf(z') 



From the above one trivially finds. 



k = k 



(6) 

(7) 
(8) 



In addition the power spectrum is sensitive to the volume 



° This expression is correct as long as the distribution of photo- 
metric errors is Gaussian, as we assume throughout this paper. 



element. Thus we must re-scale P"^°'^ by the differential vol- 
ume element with respect to the reference cosmology : c?lC\\ ■ 
Lastly, following (Tegmark 1997) and (Seo et al. 2003) we 
construct the ior each radial bin i as. 



dkt 
"8^ 



■ J dfi Gov J {k, y) {pf'' {k,n,Zi) 



}_pmod,7 ~ 



(9) 



where Cov^g is defined for every bin i according to, 



Coy J {k,ii) 



x{i) 



d^r 



n(r) 



l + n{r)P°>"'{k,fj,,Zi) 



(10) 

This is where the covariance of the power spectra is ac- 
counted for, which we assume to be diagonal in k. It has 
contributions from both sample variance and shot noise. In 
Eqs. (|9I10|) P"*"* is the measured spectra in the chosen ref- 
erence cosmology which we take as our fiducial cosmological 
model introduced in Sec. (|2.5|l . 

For the spectroscopic survey we assume that bins are 
uncorrelated. Thus the total is given by, 



2 

X3D 



E^3d(^), (11) 

i 

where the sum runs over all the bins considered. 

2.3 Angular (2D) power spectrum 

In our 2D analysis we consider the exact computation of 
the angular power spectrum of projected overdensities in a 
radial shell, 

cf = lfdk k^Poik) (^m + *r(fc))' (12) 



where 

¥t{k)= J dz cl>i(z)b{z)Diz)je(kr{z)) 
is the kernel function in real space and 

*r(fc) = 1'^' H^)f{z)D{z) 



(13) 



21^ + 21-1 



{21- 1){2£ + 1) 
{£ + !){£ + 2) 



{2l + 3){2£-l) 
ji-2{kr) 



je{kr) 



2l+l){2l + 3) 



ji+2{kr) 



(14) 



should be added to "i/l if we also include the linear Kaiser 
effect (Fisher, Scharf & Lahav 1994; Padmanabhan et al. 
2007). In turn, photo-z effects are included through the ra- 
dial selection function 0(z), see below. This model then has 
the same assumptions as the 3D spectrum from Eq. Q. 

Notice that in Eq. (|12p we are only considering den- 
sity and redshift space distortions terms. We are neglecting 
General Relativity (GR) effects as well as velocity and lens- 
ing terms, which are in our cases subdominant to the ones 
considered. Nonetheless the framework of angular auto and 
cross correlations could easily include these effects when re- 
quired (Bonvin C., Durrer R. 2011; Challinor & Lewis 2011). 
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There are A^^ angular power spectrum, one per radial 
bin. But if we want to study all the clustering information 
we should add to our observables the Nz{Nz — l)/2 cross- 
correlations between different redshift bins. These are given 
by 

^Ijdk k'p{k) (v&Kfc) + + 

(15) 

Therefore, we are considering Nz{Nz + l)/2 observable an- 
gular power spectra when reconstructing clustering informa- 
tion from tomography using A^^ bins. 

2. 3. 1 Radial selection functions 



The radial selection functions 4>i in Eqs. (|12I15|I are the prob- 
ability to include a galaxy in the given redshift bin. There- 
fore, they are the product of the galaxy redshift distribution 
and a window function that depends on selection character- 
istics (e.g binning strategy). 



(16) 



where dNg/dz is given by Eq. ([T|. We consider two different 
W{z) depending on the kind of redshifts estimation. In an 
spectroscopic redshift survey W{z) is a top hat function with 
the dimensions of the redshift bin. On the other hand, if we 
include the effect of photo-z then 

W^iz) = J dzpP{z\zp)W^{zp), (17) 

where Zp is the photometric redshift and P{z\zp) is the prob- 
ability of the true redshift to be z if the photometric estimate 
is Zp. For the photometric surveys we assume a top-hat se- 
lection W{zp) in photometric redshift and that P{z\zp) is 
gaussian with standard deviation cr^. This leads to, 



dz 



erf 



erf 



(18) 



where Zp,min and Zp^max are the (photometric) limits of each 
redshift bin considered. In the equation above and through- 
out this paper we assume cy^ is constant in redshift. 



2.3.2 Covariance matrix of angular power spectra 

The covariance between angular spectra of redshift bins ij 
and redshift bins pq is given by 



Gov 



obs ,ip y~iobs ,j q 



obs,iq ^obs,jp 



«,(ii)(P9) 



Nil) 



(19) 



where N{£) = {2£ + l)A£fsky is the number of transverse 
modes at a given £ and A£ is typically chosen to make 
Gov block-diagonal (Gabre et al. 2007; Grocce, Gabre, & 
Gaztafiaga 2011). For simplicity we consider an ideal full 
sky survey and use A£ = 1 and f^ky = 1. In this way we 
avoid correlations between different modes in the covariance 
matrix, which is diagonal with respect to £ (which is con- 
sistent with assuming the 3D covariance is also diagonal in 
k). 

Therefore, for each £ we define a matrix with N{N+1) /2 



elements, where A'^ is the number of observables discussed in 
Sec. 12.31 to account for the covariances and cross-covariances 
of auto and cross-correlations. In order to include observa- 
tional noise we add to the auto-correlations in Eq. (|19p a 
shot noise term 



obs,ij 



cy + s. 



An 



(20) 



that depends on the number of galaxies per unit solid an- 
gle included in each radial bin. We define the X2d assuming 
the observed power spectrum Ci*" correspond to our fidu- 
cial cosmological model discussed in Sec. (|2.5p . while we call 
C™°'* the one corresponding to the cosmology being tested. 



2 

X2D 



y^mod 



Gov, ^ I Ce^" 



Cmod 
t 



(21) 



Notice that each term in the sum is the product of Nz [Nz + 
l)/2-dimensional vectors C^"* where (ij) label all possible 
correlations of Nz redshift bins, and a Nz{Nz -I- l)/2 x 
Nz{Nz + l)/2 matrix corresponding to their (inverse) co- 
variance. 

Recall that we use the exact calculation of Ce using 
CAMB_sources, rather than the well-known Limber approxi- 
mation (Limber 1954). 

2.4 Nonlinear Scales 

Both X3D and xao depend sensibly on the maximum kmax 
(or minimum scale) allowed in the analysis. In this paper, 
we chose to fix kmax for all the bins and relate it to angu- 
lar scales through £max = k,naxr{z), where z is the mean 
redshift of the survey. In our fiducial cosmology we find 
r{z) = 1471 /i"^ Mpc in the redshift range 0.45 < z < 0.65 
and r{z) = 2219 /i"^ Mpc when 0.4 < z < 1.4. In addition, 
we do not consider a dependence of Imax with redshift (i.e. 
same £max for all redshift bins and their cross-correlation). 



For the largest scale we use fc„ 



10"^/iMpc- 



the 3D analysis and £m.in = 2 in the angular case. We have 
not found any significant dependence on Kmin or £min ■ 



2.5 Cosmological model and growth history 

We assume the underlying cosmological model to be a flat 
AGDM universe with cosmological parameters w = —1, h = 
0.73, Us = 0.95, Qm = 0.24, ^i, = 0.042 and ag = 0.755. 
These parameters specify the cosmic history as well as the 
linear spectrum of fluctuations Pq. In turn, the growth rate 
can be well approximated by, 



(22) 



and 7 = 0.545 for AGDM. Gonsistently with this we obtain 
the growth history as 



D{z) 



exp 



f{z) 
l + z^ 



dz 



(23) 



(where D is normalized to unity today). The parameter 7 is 
usually employed as an effective way of characterizing mod- 
ifled gravity models that share the same cosmic history as 
GR but different growth history (Linder 2005). In part of 
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our analysis we focus in ACDM models and assume the GR 
value 7 = 0.545. We deviate from this in Sec. 13.1.21 were we 
take 7 as a free parameter independent of redshift. 

2.6 Likelihood analysis 

In order to find constraints on cosmological models we in- 
tegrate over the space of parameters defining the model, 
finding the value of the likelihood given by 

-21og£(Xx', (24) 

where we approximate the likelihood as Gaussian in the 
power spectra. Given the prior ■& on the parameters one 
defines a probability for each sampled point i in parameter 
space given by 

V{i) oc C(i) X (25) 

Finally, the mean and covariance matrix of the parameters 
is obtained from 

Pa = ^^(i)Pa(i) (26) 

i 

^(Pa,P,) = J2^(^)(Mi)-Pa){Pb{i)-pt), (27) 

i 

where Pa{i) is the value of the parameter a in the grid point i, 
Pa is the mean value and E(p^_pj^) is the covariance between 
parameter a and b. In Eqs. (|26I27|) ^(i) is normalized to 
unity over the grid. In addition we assume flat priors. 

By construction the likelihood peaks at the fiducial 
value considered in the analysis. In all our studies we have 
chosen wide prior limits and therefore have found no de- 
pendance with these limits, and find the mean agrees with 
the fiducial value and the posteriors are quite Gaussian. 
Then in the case of only one nuisance parameter p, solving 
X'^(p) — 1 = gives the same variance as likelihood sampling 
which allows us to speed up constraints considerably. 



where 5* is the subspace of parameters we are interested in. 
If this subspace correspond to only one parameter, then the 
FoM is the inverse of the square root of the variance of the 
corresponding parameter. Thus we have the following cases, 

• FoMo„: Constraints on i}m,, with other parameters 
fixed at fiducial values. 

• FoMi, and FoM-yi bias and 7 constraints when marginal- 
ized over 7 and bias, respectively. Other parameters are fixed 
at their fiducial values. 

• FoMj,.y: Joint constraint on bias and 7, with other pa- 
rameters fixed at fiducial values. 



3 RESULTS 

In this section we present the forecasts on flm {bias fixed) 
and b and 7 {bias free) from the measurement of either spa- 
tial or angular power spectra in the spectroscopic survey 
described in Sec. 12.1.11 Next we perform the bias fixed anal- 
ysis in the narrow-band photometric survey with accurate 
photo-z discussed in Sec. 12.1.21 and the broad-band photo- 
metric survey defined in Sec. 12.1. SI Notice that despite pho- 
tometric redshift errors large scale redshift space distortions 
can be measured in photometric surveys for binned data 
(Nock, Percival, & Ross 2010; Crocce, Cabre, & Gaztafiaga 
2011; Crocce et al. 2011), albeit with possible large error 
bars. Nonetheless for photometric surveys we concentrate 
the bias fixed case only. 

All the analyses introduced above have been done for 
three different kmax = {0.05, 0.1, 0.2} /i Mpc"^ (with cor- 
responding imax as detailed in Sec. 12. 4p . and several bin 
configurations (see Table [2] and [3]). We then study for which 
redshift bin width the information obtained using angular 
power spectra (quantified by the FoM of Sec. 12. 7p are simi- 
lar to those derived from the 3D power spectra. 



2.7 Figures of Merit 

We consider two different analyses in order to compare 3D 
clustering with 2D tomography including all the auto and 
cross-correlations between redshift bins. 

On the one hand, a bias fixed case, in which we only 
vary Sim (which affects both the shape and the amplitude 
of the power spectrum, and can be constrained as if we had 
a good knowledge of the bias prior to the analysis). 

On the other hand we consider a bias free case, in which 
only b and 7 (hence / through Eq. ((22])) are allowed to vary. 
This changes the (anisotropic) amplitude of the power spec- 
trum, but not the underlying shape. This case is virtually 
the same as the standard analysis of redshift space distor- 
tions (White et al. 2009; Ross et al. 2011). For this case 
we had to adapt CAMB_sources slightly, see the discussion 
in Appendix A. 

To make the comparison quantitative we define a figure 
of merit (FoM) based on the covariance matrix E, 



3.1 Spectroscopic redshifts 

3.1.1 Bias fixed case 

Top panels of Fig. [2] show the FoM on Qm for different kmax 
and imax = r{z)kmax as a function of the number of red- 
shift bins N^, in which we divide the full survey volume (see 
Table [2]). Here dashed lines are results from fitting the 3D 
power spectrum according to Eqs. (|9I11[I . while solid are 
from the 2D tomography including all the auto and cross 
correlations of bins, as in Eq. (|2H) . The left (right) panel 
corresponds to the low (high) shot noise case as defined in 
Tabled] 

As expected we find that the FoM increases for increas- 
ing kmax, imax. lucludiug morc modes to the adds more 
information to our analysis and therefore results in better 
constraints. We also see that FoMn„ from the 3D analysis 
only show a marginal dependance on the bin configuration. 
This is because the per redshift bin is roughly propor- 
tional to the volume of the redshift shell, see Eq. (|9]l . Thus, 
increasing the number of bins at the expense of decreasing 
their volume keeps the FoMo„ unchanged. We obtained the 
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Figure 2. Spectroscopic survey & bias fixed. Top panels show FoMf2^(2D) and FoMfj^ (3D) as a function of the number of bins in 
which we divide the survey for the analysis (left panel for a low shot- noise survey and right to a high shot noise). Dashed line corresponds 
to the 3D analysis, dotted to the 2D tomography using only auto-correlations and solid to auto plus cross correlations. Different colors 
correspond to different minimum scales, as detailed in the bottom panel inset labels. Bottom panels show the ratio of FoMq^ (2D) (auto 
plus cross) and FoMfj^ (3D) as a function of the bin width Ar normalized by the minimum scale assumed in the 3D analysis. Remarkably 
the recovered constraints from full tomography match the 3D ones for Ar ~ for all . We note that different lines in the bottom 

panels are truncated differently merely because we have done the three kmax cases down to the same minimum Ar. 



same result for all the cases studied in this paper, as long as 
Pg does not change abruptly with redshift. Thus from now 
on we will only refer to the 3D results in the whole survey. 

This picture changes for the 2D tomography. Here the 
transverse information is fixed once imax is set {2£+l modes 
per £ value up to Imax)- As we increase the number of nar- 
rower bins Nz (with fixed total redshift range) we have sev- 
eral effects: 

(i) Decreasing the number of galaxies per bin increases 
the shot noise per bin 

(ii) Increasing the number of bins so that they are 
thinner proportionally increases the signal auto power 
spectrum in each bin (there is less signal power suppression 
due to averaging along the radial direction). 

(iii) When we split a wide redshift bin in two, we double 



the number of angular auto power spectra (transverse 
modes). This results in a larger FoM because the signal 
to noise in each bin remains nearly constant (the shot 
noise and signal in each bin both increase proportionately). 
This gain is illustrated by the dotted line in Fig. [21 
which corresponds to the FoM produced by just using 
auto-correlations. For even narrower redshift bins the bins 
will become correlated and the gain will saturate, but this 
is not yet the case in our results as the redshift bins are still 
large compared to the clustering correlation length. In the 
limit in which all modes of interest are very small compared 
to shell thickness and they are statistically equivalent, for 
a single power spectrum amplitude parameter one expects 
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Figure 3. Spectroscopic survey & bias free. Top panels show the combined b — 7 constraint resulting from 3D clustering (dashed lines) or 
2D tomography considering as observables only auto correlations in redshift bins (dotted lines), or adding to this the cross-correlations 
(solid lines). The x-axis corresponds to the number of radial bins considered in the analysis. Different colors label different minimum 
scales assumed (same values and labels as in Fig. [2]| . Middle and bottom correspond to individual b or 7 constraints after marginalization 
over 7 or b respectively. As for the bias fixed we find that 3D information can be recovered but now the role of radial modes is much for 
important because RSD (our bias free case) relies on the relative clustering amplitude of radial and transverse mode. 



FoM= 1/0- oc ^/]V^, as obtained in Fig. [2] for low iV^ 0. 

(iv) When we increase the number of narrower bins, we 
also include information of radial modes by adding the 
cross-correlation between different redshift bins (illustrated 
by the solid line in Fig. [2] that corresponds to the total 
FoM from auto plus cross-correlations). Note how adding 
the cross-correlations to the autocorrelations (solid lines in 
Fig. [2]) only increases the FoM moderately as compared to 
the autocorrelation result (dotted line). This reflects the 
fact that there are fewer radial modes than transverse ones, 
while much of the Qrn constraint comes from the shape of 



^ A similar effect can be seen on Figs. 8 and 9 of (Ross et al. 2011) 
in the context of RSD constraints in a broad band photometric 
survey. In their Fig. 8 the constraint in / <Tg saturates when they 
consider only one redshift bin. However the error on /erg from the 
combined measurements on several bins does not saturate (Fig 9) . 



P{k) that is isotropic. 

(v) As shown in Fig. [H the 2D FoM can exceed the 3D 
FoM. This happens because the 3D analysis is limited by 
construction to a maximum number of modes, given by 
kmax, while in 2D we only limit the analysis to I < Imax 
and we can formally exceed the maximum number of nar- 
row redshift bins, as explained in point (ii) and (iii) above. 
But in reality, these additional modes are not necesarilly in- 
dependent and they could well be in the non-linear regime, 
so it is not clear to what extent we can use them to increase 
the FoM. As we want to restrict our analysis to A: < kmax 
we should not use redshift bins that are smaller than A^„. 

The bottom panels of Fig. [2] show the ratio of the 2D 
and 3D FoM's against the bin width (instead of N^), now 
normalized by the minimum scale used in the 3D analysis 
^min = k'^^ (for three different kmax as before). We find 
FoM(2D) ~ FoM(3D) when A5tS„ ~ Ar for all More 



© 2012 RAS, MNRAS 000, 1-10 



9 

Low Shot Noise High Shot Noise 




Figure 4. Spectroscopic survey & bias free. Top panels show the ratio between combined FoM;,-^ (2D) (auto plus cross correlations) and 
FoMi,-y (3D) with respect to A^^^^ = 2n/kmax, normalized by the mean width of the redshift bins Ar in the analysis. Middle and bottom 
panels show the same but for ratios of FoMj, and FoM-y, respectively. We show results for Low Shot Noise and High Shot Noise in left 
and right panels, respectively. To reconstruct RSD information in practice, one need bins slightly smaller than A^^^^^. 



precisely: 

Ar^cAz/H{z)~O.SXTi„. (29) 

Basically this means that the 3D clustering information 
is recovered once the binning is such that the radial bin 
width equals the minimum scale probed in the 3D analysis. 
In this case one is able to constrain the parameters without 
loss of information compared to a three dimensional analy- 
sis, though the actual range of scales around kmax that are 
used in the 2D analysis may be slightly different from the 
ones used in the 3D analysis. 

Note that as mentioned in point (v) above, we can only 
really trust our results for the 2D FoM up to the limit in 
which they are equal or smaller than the 3D FoM, i.e., in the 
range in which the width of redshift bins is greater or similar 
than A^„. To use the smaller scales we first need to explore 

© 2012 RAS, MNRAS 000, 1-10 



to what extent we can model the non-linear 2D clustering to 
improve the FoM. We are currently investigating this issue 
(Asorey et al., in prep.) 

Lastly, note that including shot noise does degrade the 
FoM as shown in the right panel of Fig. [2] However this does 
not change the conclusions above. 

3.1.2 Bias free case 

We now turn to the bias free case where we assume we know 
perfectly the shape of the power spectrum so that all the 
parameters are fixed at their fiducial values except the bias 
b and the growth index 7. 

In Fig. [3] we plot the combined FoM obtained for bias b 
and growth index 7, and the FoM of each of these 2 parame- 
ters marginalized over the other, as a function of the number 
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Figure 5. Narrow band photometric survey (PAU-like) & bias fixed. Top panels show figures of merit FoMn^ (2D) (auto plus cross 
correlations) and FoMq^ (3D) with respect to the number of bins for kmax = {0.05,0.1,0.2} hMpc^^ (red, violet and orange colours). 
We plot 2D figures of merit with solid lines and 3D figures of merit using dashed lines. Bottom panels show the ratio between both 
figures of merit with respect to minimum scale used in 3D analysis, = 2TT/kmax, divided by the mean width Ar of the redshift 

bin. We conclude that we get similar constraints from 2D and 3D analysis when Ar is close to A^^^ji ^-i^d that in terms of bind width 
optimization an spectroscopic and photometric analysis are almost identical. 



of redshift bins considered in the analysis (for a fixed survey 
redshift range 0.45 < z < 0.65). As in Fig. [S] dashed line 
corresponds to the 3D analysis, dotted hne to the 2D tomog- 
raphy using only auto-correlationj3 B^nd solid line to the full 
2D case where we add auto and cross angular correlations. 

We find a similar trend for the evolution of the dif- 
ferent FoM of the 7 and 6 parameters (either combined or 
marginalized) than when varying $1™,. Constraints given by 
spatial power spectrum are stable, while constraints from 
projected power spectrum in the bins increases with the 
number of bins in which we divide the survey. However there 
is a substantial difference in regards to the contribution of 
radial modes. Now the contribution of cross-correlations is 
very large (compare solid to dotted lines in the left panel of 



" We note that we refer here to observables. The covariancc of 
the auto-correlations does include cross-correlations of redshift 
bins, see Eq. I I19I1 . 



Fig. ll9|) . Iir fact, without cross-correlations we do not recover 
all the 3D information. This is because redshift space distor- 
tion information (i.e. our bias free case) is based in the rela- 
tive clustering amplitude of modes parallel and transverse to 
the line of sight. The contribution from radial modes is much 
more evident for the 7 constraint (FoM^, and then FoMi,.y) 
because 7 is basically what quantifies this relative clustering 
amplitude (in addition / = i}{z)'^ depends on redshift while 
we assume bias does not). 

As we have done with FoMn,„ we show in Fig. |4] the 
dependence of the ratios between 2D and 3D FoM with 
respect to A^„/Ar. We find that both analyses produce 
the same constraints when the mean redshift bin width is 
slightly smaller than (and we use auto and cross 2D 

correlations). Comparing these results with the bias fixed 
case, it seems that for the RSD probe we need to extract 
more radial information. In this case: 

Ar = cAz/H{z) ~ 0.6 A^^„ (30) 
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Figure 6. Broad band photometric survey (DES-like) & bias fixed. Top right panel shows the figures of merit FoMq^ (2D) and FoMq^ 
(3D) with respect to the number of bins Nz used in the analysis, for kmax = {0.05, 0.1, 0.2} /i Mpc~^ (red, violet and orange colours, 
respectively). 2D FoM are plotted with solid lines and 3D with dashed lines, and we only consider low shot noise. Bottom left panel show 
the ratio of both figures of merit with respect to A^^^ = 2-K/kmax divided by Ar. The equivalence of the recovered FoM now changes 
for different kmax- However when this ratio is plotted with respect to the comoving scale of photo-z, cr^ (normalized by Ar) the different 
Amin lines cross each other for Ar ~ 2ar ■ This implies that is the relative values of Ar and Or what sets the equivalence of 3D and 2D 
tomography. In particular, for a DES-likc survey one recovers the 3D constraints from 2D analysis using 5 rcdshift bins. 



as compared to 0.8 in Eq. (|29|) . This means that we have 
to include more radial bins when developing the fit to an- 
gular correlations than when only fitting $1^ if we want to 
match the constraints from 3D clustering. This in practice 
corresponds to using slightly narrower redshift bins. This 
may also result in more information being included from ra- 
dial modes with k > kmax, though a detailed analysis of the 
implications of this is beyond the scope of the current paper. 

3.2 Photometric redshifts 

In this section we show how the results found in the pre- 
vious section extend to the photometric surveys detailed in 
Sec. I2.1.^ and [?.1.3I For concreteness we will only consider 
the bias fixed study where all cosmological parameters are 
fixed at their fiducial values except for flm- 

3.2.1 Narrow-hand photometric survey (PAU-like) 

In top panels of Fig. [S] we show the Q.m constraints {bias 
fixed case) from 3D and 2D analysis (dashed and solid lines 
respectively) in a narrow band photometric survey with 
(Jz = 0.004. In bottom panels we show how the ratio be- 
tween 2D and 3D FoM depends on the ratio between the 
minimum scale of the 3D analysis and the mean comoving 
width of radial shells. 

We find basically the same result as in the spectroscopic 



survey. Constraints from a projected or unprojected analy- 
sis are equivalent when the mean width of the radial shells 
(set by our binning strategy) is equal to the minimum scale 
considered in 3D analysis A^„. The absolute value of each 
FoM is degraded with respect the FoM reached with an spec- 
troscopic survey because photo-z errors dilute clustering in 
the radial direction. This broadens the selection functions 
in the 2D analysis and introduces a cut off already at quasi- 
linear scales in the 3D P{k). In both cases the consequence 
is that signal to noise reduces and thus errors of observables 
degrade. But if we compare Fig. [2] and Fig. [5] we see that 
the spectroscopic survey and a photometric one with very 
accurate redshifts are almost indistinguishable in terms of 
bin width optimization. 



3.2.2 Broad-band photometric survey (DES-like) 

We now consider a deep survey (ias < 24) with redshifts 
estimated by photometry with broadband filters (cr^ = 0.1), 
and use the full catalogue with 0.4 < 2; < 1.4. We obtain the 
FoM for Qm shown in the top left panel of Fig. [S] 

Now the large photo-z error removes most of the radial 
information, thus all FoMn^^ are degraded with respect to 
spectroscopic and narrow-band photometric surveys. In ad- 
dition, we find that FoMn^ saturates with the number of 
redshift bins included in the survey for 
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feet is produced by the overlapping between true galaxy dis- 
tributions at different bins induced by photo-z transitions. 

We also find that the configuration in which spatial and 
projected analysis constrain Qrn equally corresponds to the 
same number of bins for all the kmax considered. Therefore, 
as we can see in bottom left panel of Fig.[Hl the scale given by 
^min is not ruling the dependencies. Instead it is the scale of 
the photometric redshifts which is affecting both clustering 
analyses. This is shown in the right panel of Fig. [S] where 
we plot the ratio of figures of merit (2D vs. 3D) against a 
new scaling : ar/Ar. We find that for a DES-like case, with 
the assumption of cr^ = 0.1, one needs roughly 5 bins for 
the 2D tomography to optimally recover the 3D clustering 
information. This corresponds to: 

Az ~ 2a^. (31) 

With a lower the number of bins will increase. 



4 CONCLUSIONS 

In this paper we have studied the redshift bin width that 
allows us to recover the full 3D clustering constraints from 
tomography of angular clustering (i.e. the combination of 
all the auto and cross correlations of redshift bins). We ex- 
plore three surveys with different properties: a spectroscopic 
and a narrow band photometric survey in a redshift range 
0.45 < z < 0.65, and a deeper broadband photometric sur- 
vey that covers redshifts in the range 0.4 < z < 1.4. We 
have considered how well we can recover the shape of the 
power spectrum by allowing $1™ to be free and fixing the 
amplitude of clustering, including bias. We call this the bias 
fixed case. We have also explored how to recover the infor- 
mation from redshift space distortions (RSD), by measuring 
the anisotropic amplitude of the power spectrum allowing for 
both a free bias and a free growth index. This is the bias free 
case. We restrict our study to quasi-linear scales and we only 
consider scales above some minimum scale — 2Tv/kmax, 
where k < kmax and kmax is either 0.05, 0.1 or 0.2 /iMpc^^. 
In angular space this corresponds to / < Imax — kmaxr{z), 
where r{z) is the radial distance to the mean redshift bin. 

The 3D analysis has almost no dependance on the num- 
ber of redshift bins because radial modes are already in- 
cluded in each bin. In contrast the 2D tomographic analysis 
depends strongly on the number of bins (or equivalently on 
redshift bin widths), since broad bins average down trans- 
verse power on scales smaller than the bin width, and it 
is only by using multiple thin shells that radial modes are 
included. 

For the bias fixed case in the spectroscopic survey we 
have found that we recover all the information with 2D 
tomography when the width of the redshift bins that we 
use to do the tomography is similar to the minimum scale 
used in the 3D observables, A^„ . More precisely we find 
that the optimal bin width is (see Fig. [2] and Eq. (|29|) '): 
Ar = cAz/H{z) ~ 0.8 A^„. In addition most of the 2D 
constraints come from autocorrelations. 

When studying RSD, i.e. in the bias free case, we see 
that if we want to recover the 3D constraints we need ra- 
dial shells which are slightly smaller, i.e. Ar ~ 0.6A^„ (see 



Fig- Bli which means that we would need more bins than 
in the case in which we just want to measure the shape of 
P{k). In addition we find necessary to include in the ob- 
servables the cross correlation between redshift bins. This 
is expected because in the RSD case we are comparing the 
clustering in radial and transverse direction to the light of 
sight: information from radial modes should be more impor- 
tant than in the case in which we just study information in 
the isotropic shape of the power spectrum. Also note how 
we can not recover the 3D information from RSD when we 
just use autocorrelations (see dotted line in Fig. |3}. 

We found that in the bias fixed case, the narrow-band 
photometric survey is almost equivalent to an spectroscopic 
survey, and we therefore reach the same conclusions with re- 
spect to the optimal bin width for the tomography of galaxy 
counts. In the case of a deeper broadband photometric sur- 
vey we find that the typical uncertainty in photometric red- 
shifts (Tz severely limits the accuracy of the radial informa- 
tion for both 3D and 2D cases. In this case the information 
recovery does not depend strongly on A^„, because this is 
smaller than the scale corresponding to the photometric red- 
shift accuracy, i.e. ca^jHiz) > A^„. The optimal redshift 
bin width in this case is simply given by Az ~ 2(T2. 

For a redshift range 0.4 < z < 1.4 and = 0.1 (DES- 
like survey) we find that we will need only 5 redshift bins to 
constrain f2m using tomography with the similar precision 
than a full 3D analysis of the survey. In comparison, for a 
PAU-like survey with — 0.004 and kmax = 0.1 we need 
about 44 redshift bins of width Az ~ 0.023 each. 

We conclude from our analysis that it seems possible 
to recover the full 3D clustering information, including RSD 
information, from 2D tomography. This has the disadvan- 
tage of needing a potentially large number of redshift bins, 
and correspondingly large covariance matrices between ob- 
servables. But it has the great advantage of simplifying the 
combination with WL and of just using observed quanti- 
ties, i.e. angles and redshifts, avoiding the use of a fiducial 
cosmology to convert angles and redshifts into 3D comov- 
ing coordinates. In practice, probably both types of analysis 
should be used to seek for consistency. 
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and cross-correlation is given by: 



Cf = hicf + 26./,C; + ffcf (Al) 



CI- 



(A2) 



where bi is the bias of the bin i and fi is the growth rate 
given by Eq. (|22[) . evaluated at the mean redshift of the bin i. 
This factorization assumes f{z) does not vary much within 
the redshift range of the bin (neither b). We have tested 
this assumption using the exact CAMB_sources evaluation or 
the reconstruction of Eqs. HA1IA2|I and found an excellent 
match for the bin widths considered in this paper. 

Using the observed Ce and solving a linear set of equa- 
tions using different values for bi we can store the value 



of a 



a (2) ^ij (2) ^ij {2') ^ii (4) 



c;' <°)and c;' are obtaining by excluding RSD in Ci. 
Then, we sample b and 7 space using these factors and the 
reconstruction given by Eqs. HA1|) and (|A2|I obtaining C™"'' 
in parameter space. 

In the reconstruction we assume the underlying value 
of Sim = 0.24 given by our reference cosmology while the 
growth factor D{z) is included in the integrals that are con- 
tained in the cosmic history dependent factors C^"' . 

This paper has been typeset from a TlgX/ I^T^eX file prepared 
by the author. 



-lij (0) 



and C, 



ij (4) 



The values of 



APPENDIX A: MODIFYING CAMB_SOURCES TO 
SAMPLE GROWTH RATE AND BIAS 

In order to consider the bias free case we had to modify 
CAMB_sources to accept as (independent) inputs bias and 
growth rate (parameterized through 7 as in Eq. (|22p ). In 
addition this case does not involve changes in the shape of 
the real space spectrum, thus one should be able to sample 
parameter space without the need to compute the transfer 
functions at each point of parameter space. 

To fulfil these needs we have factorized the terms in 
our observables that depend on the cosmic history (for our 
reference cosmology) from those that depend on the bias b 
and growth index 7. The factorization in the case of auto 
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